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ABSTRACT 

High  elevation , panoramic  infrared  photography  taken 
from  a U-2  aircraft  was  evaluated  to  determine  its  effec- 
tiveness in  quantifying  annual  mortality  of  lodgepole  pine 
caused  by  the  mountain  pine  bettle  in  a 520, 650-acre  outbreak 
in  portions  of  the  Beaverhead  and  Gallatin  National  Forests 
in  Montana.  A multistage,  variable  probability  design  using 
probability  proportional  to  size  (PPS)  at  three  levels  was  used 
throughout  the  survey. 

The  results  suggest  that  panoramic  photography  can  be 
effectively  used  to  provide  precise  estimates  of  mountain  pine 
beetle-caused  annual  mortality  of  lodgepole  pine  over  large 
areas.  Mortality  of  lodgepole  pine  in  1977  was  estimated  at 
1,891,510  ±195,805  trees,  and  27,001  M ±3682  M cubic  feet  of 
volume.  These  estimates  differed  somewhat  from  a multistage 
aerial  photo  survey  also  conducted  in  1978,  the  difference 
ascribed  to  different  areas  of  coverage.  Physical  characteris- 
tics of  this  unconventional  format  are  described,  and 
suggestions  for  its  future  use  in  similar  bark  beetle  damage 
surveys  are  recommended. 
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1.  Portions  of  the  Beaverhead  and  Gallatin  National 
Forests  showing  the  location  of  the  OB  frames  (±  UO°)  and 
distribution  of  the  mountain  pine  beetle  (in  red)  as 
determined  during  a 1978  aerial  sketchmap  survey. 
Nadir  is  shown  as  a circle  to  the  left  of  the  frame  numbers. 


INTRODUCTION 

Sound  data  on  tree  damage  and  mortality 
caused  by  destructive  forest  insects  and 
diseases  is  required  for  pest  management 
decision  making  at  all  planning  levels.  During 
the  past  two  years,  Forest  Insect  and  Disease 
Management  has  conducted  multistage  sur- 
veys to  estimate  annual  mortality  of  lodgepole 
( Pinus  contorta  var.  latifolia  Engelm.)  and  ^ 
ponderosa  pine  ( Pinus  ponderosa  Laws)y 
caused  by  the  mountain  pine  beetle  (jpendroc£‘ 
tonus  ponderosal  Hopkins).  These  surveys, 
conducted  in  portions  of  the  Targhee  National 
Forest,  Idaho,  and  the  Black  Hills  National 
Forest,  South  Dakota  and  Wyoming,  have 
provided  precise  estimates  of  mountain  pine 
bettle-caused  tree  mortality  over  large  areas 
with  reasonable  sampling  errors  (Klein  et  al. 
1979,  Hostetler  and  Young  1979a,b).  Flexi- 
bility of  the  general  survey  design  allows  for 
application  to  small  areas  as  well  as  large 
areas,  e.g.,  a National  Forest. 

Medium-scale,  reconnaissance  photography 
in  a panoramic  format  taken  from  high 
altitudes  (65,000  feet)  has  recently  been  made 
available  for  resource  management  applica- 
tions. This  photography  has  an  advantage  over 
conventional  aerial  photography  in  that  large 
areas  can  be  photographed  in  a relatively  short 
period  of  time.  The  panoramic  film  can  then  be 
used  for  preliminary  classification  as  well  as 
direct  sampling,  thus  eliminating  need  for 
aerial  sketchmap  surveys.  The  photography  is 
acquired  by  the  Itek  KA80A  optical  bar  (OB) 
camera,  and  has  been  used  to  estimate  volumes 
of  dead  timber  on  the  Clearwater  National 
Forest,  Idaho  (Duggan  et  al.  1977),  and  to 
stratify  and  quantify  insect-killed  trees  in 
California  (Klein  et.  al.  1978). 


DESCRIPTION  OF 
TEST  SITE 

The  test  site  was  in  portions  of  thejGallatin 
and  Beaverhead  National  ForestjQm  south 
centrj^^ontanajTfigure  T).  The  area  is 
com  jguseH  of  parts  of  the  Madison  and  Gallatin 
Mountain  ranges,  separated  by  the  Gallatin 
River,  and  lying  between  the  Yellowstone  and 
Madison  Rivers.  Elevations  range  from  slight- 
ly over  5000  feet  in  the  Gallatin  River  Valley  to 
more  than  11,000  feet  in  the  Madison  Range. 
The  generally  steep  terrain  consists  of  U- 
shaped  glacial  valleys,  sharp  ridges,  almost 
vertical  headwalls,  and  treeless  peaks.  The 
dominant  conifer  is  lodgepole  pine,wifch-  lesser 
amounts  of  Douglas-fir  [Pseudat&tfga  menziesii 
var.  glauca  (Beissn.)  Franco],  Engelmann 
spruce  (Picea  engelmannii  Parry  ex  Engelm.), 
subalpine  fir  [Abies  lasiocarpa  (Hook.)  Nutt.], 
whitebark-limber  pine  ( Pinus  spp.),  and  pon- 
derosa pine. 

The  present  mountain  pine  beetle  outbreak 
began  in  the  late  1960’s.  These  infestations 
occur  in  portions  of  the  Gallatin  and  Hebgen 
Lake  Ranger  Districts  of  the  Gallatin  National 
Forest,  and  in  the  Jack  Creek  and  West  Fork 
drainages  of  the  Madison  River  on  the  Beaver- 
head National  Forest.  During  1978,  upwards 
of  60  trees  per  acre  were  killed  in  the  Hebgen 
Lake  Ranger  District  (McGregor  et  al.  1978). 
In  many  locales  the  beetle  has  moved  into 
higher  elevations  and  is  attacking  whitebark 
and  limber  pines. 


OBJECTIVE 

Objective  of  this  survey  was  to  evaluate  the 
capabilities  of  KA80A  OB  color  infrared 
photography  for  measuring  annual  mortality 
of  lodgepole  pine  caused  by  the  mountain  pine 
beetle.  In  addition,  an  attempt  was  made  to 
produce  an  accurate  map  showing  the  distribu- 
tion and  relative  intensity  of  tree  mortality. 
These  maps  could  be  integrated  with  a geo- 
graphic data  base  and  serve  as  a management 
tool. 
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CHARACTERISTICS 
OF  PANORAMIC 
PHOTOGRAPHY 

A single  panoramic  frame,  produced  by  a 
120°  sweep  of  the  KA80A  Itek  camera, 
measures  4.5  x 50  inches  and  covers  an  area  of 
approximately  2.6  miles  (intrack  at  nadir)  by 
43  miles  (crosstrack).  The  scale  at  nadir  is 
approximately  1:30, 0002  assuming  average 
terrain  elevations  of  5000  feet.  The  scale 
decreases  outward  from  nadir  and  reaches  a 
minimum  of  1:60,000,  because  of  the  pan- 
oramic angle.  Due  to  the  decreasing  scale, 
increasing  image  obliquity,  and  color  attenua- 
tion problems  near  the  film  edge,  the  effective 
working  angle  was  judged  to  be  80°  or  ± 40° 
from  nadir.  These  characteristics  are  unique 
to  panoramic  photography,  and  though 
present  in  metric  frame  camera  systems,  they 
generally  are  not  as  severe.  Additionally, 
panoramic  photography  presents  geometric 
problems  that  result  from  panoramic  and 
topographic  distortions,  and  from  variations  in 
aircraft  attitude.  Panoramic  distortion  in  this 
survey  was  removed  by  use’of  equal-area  grid 
overlays  (nadir  scale  1:30,000).  Panoramic 
distortion  as  well  as  other  errors  can  be 
removed  by  map-registered  grid  overlays 
which  permit  relatively  accurate  transfer  of 
data  from  photo  to  maps.  The  systems  that 
provide  both  of  these  computer-generated 
grids  were  developed  by  the  Forest  Service 
WO  Geometronics  Development  Group3.  Map- 
registered  grids,  although  more  accurate, 
require  more  effort  to  produce,  and  are 
considerably  more  expensive  than  equal-area 
grids. 


2 The  scale  at  sea  level  is  1:32,500. 

3 Liston.  R.L.  1978.  Grid  generation  cost  and  pro- 

duction rate  analysis  for  the  Northern  California 

Test  Project.  Multilithed.  15  pp. 


METHODS 

Survey  Design  in  Brief 

A multistage,  variable  probability  design 
was  used  for  this  survey.  A total  of  79  alternate 
panoramic  frames  encompassing  the  infesta- 
tion were  obtained  during  the  photo  mission. 
These  were  divided  into  160-acre  squares 
using  equal-area  transparent  grids  (figure  2). 
Experienced  aerial  photo  interpreters  using 
specially  adapted  microfiche  viewers  made 
quick  estimates  of  the  number  of  faders  within 
each  160-acre  plot.  This  first  cut  produced  a 
list  from  which  120  160-acre  plots  were 
selected  for  detailed  photo  interpretation  by 
probability  proportional  to  size  (PPS).  These 
plots  were  then  further  divided  into  sixteen  10-acre 
squares  using  a second  series  of  equal-area 
grids.  All  faders4  within  each  10-acre  plot  were 
then  counted  in  stereo.  From  this  second  list, 
thirty  10-acre  plots  were  then  selected  by  PPS 
for  ground  truth.  However,  only  20  plots  were 
actually  examined  because  of  time  and  person- 
nel constraints.  Ground  crews,  viewing  the 
photography  in  stereo,  located  the  plots,  estab- 
lished plot  boundaries,  and  systematically 
recorded  all  faders. 

Aerial  Photography 

Photography  was  taken  July  25,  between 
12:12  and  12:59  pm  MDT  (47  minutes),  by  a 
U-2  aircraft  flying  from  Moffett  Field,  Califor- 
nia. The  aircraft  flew  at  19,800  M (65,000  feet) 
at  400  nm  (460  miles)  per  hour.  The  aircraft 
was  to  fly  three  parallel  flight  lines  at  an  inter- 
val of  17  nm  (19.6  miles).  High  definition  Aero- 
chrome  infrared  reconnaissance  film  (SO-131) 
was  taken  with  the  KA80A  Itek  OB  camera. 
Two  hundred  and  fifty-one  frames  were  taken 
in  the  stereo  mode.  Theoretically,  the  outbreak 
should  have  been  covered  with  79  alternate 
(nonstereo)  frames,  but  due  to  aircraft  align- 
ment problems  and  the  inherent  voids  between 
adjacent  equal-area  grids,  complete  coverage 
was  not  achieved.  The  target  area  was  cloud- 
free  and  the  quality  of  photography  was 
judged  excellent  by  NASA  standards5. 

4 A fader  is  a mountain  pine  beetle-killed  lodgepole 
pine  whose  crown  discolors  the  season  (1978)  fol- 
lowing attack  (1977).  An  (old)  fader  is  a tree 
attacked  and  killed  the  previous  year  (1976). 

5 Flight  Summary  Report  1165.  NASA  Ames 
Research  Center,  Moffett  Field,  CA,  July  25, 
1978.4  pp. 
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2.  Grid  overlays  of  1:30,000  scale  (at  nadir)  used  with  OB 
panoramic  photography. 

Top:  Equal-area  g', rid  used  in  this  survey. 

Bottom:  Map-registered  grid  used  in  the  Northern  Califor- 
nia study  (Klein  et  al.  1978).  This  flight  was  floum  on  a 
diagonal  course,  which  is  not  recommended. 


During  an  earlier  study  in  which  OB  and 
frame  camera  systems  were  compared,  the 
color  quality  of  the  OB  imagery  was  judged  to 
be  relatively  unsatisfactory  (Klein  et  al.  1978). 
On  August  8, 1978,  prior  to  the  first  film  dupli- 
cation run,  selective  color  filters  were  used  to 
enhance  the  appearance  of  the  faded  tree 
crowns  (#10  red)  and  to  subdue  the  normal 
cyan  overtone  (#20  yellow).  This  process  im- 
proved the  overall  color  balance  but  may  have 
resulted  in  a slight  degradation  of  spatial  reso- 
lution. The  overall  quality  of  the  imagery  still 
was  judged  excellent. 

Indexing  the  Photography 
The  panoramic  photography  was  obtained 
from  Moffett  Field  on  September  17,  1979. 
Prior  to  the  initial  classification  process, 
1:30,000  scale  equal-area  transparent  overlays 
were  temporarily  attached  to  the  emulsion 
(top)  side  of  those  alternate  frames  that 
covered  the  infestation  (figures  2,  3).  In  order 
to  locate  each  frame  on  the  map  (figure  1),  it 
was  necessary  to  construct  a second  transpar- 
ent grid  at  map  scale,  i.e.,  1:126,720  (1  inch  = 2 
miles).  Nadir  of  each  frame  was  then  located 
and  transferred  to  the  map.  Two  crosstrack 
points  that  could  be  identified  on  both  the  pho- 
tography and  the  map  were  then  transferred  to 
the  map.  These  conjugate  points  were  then 
located  on  the  large  grid,  transferred  to  the 
small  grid,  and  then  matched  to  their  conju- 
gate points  on  the  map.  In  the  event  the  points 
did  not  match  (due  to  elevational  displace- 
ment) they  were  equally  apportioned  on  the 
map. 


3.  Portion  of  OB  frame  near  nadir  overlaid  with  transpar- 
ent equal-area  grid  (35mm  Kodachrome  enlargement). 


Classification 

First  step  in  the  initial  classification  process 
was  to  minimize  the  actual  photointerpretation 
by  scanning  frames  over  a conventional  light 
table  and  eliminating  all  cells  not  requiring 
examination,  including  areas  of  non-forest  or 
non-host  type,  and  those  duplicated  in  the  endlap 
regions  (figures  4,  5).  Interpreters  were  then 
able  to  concentrate  only  on  those  cells  having 
lodgepole  pine  and  a likelihood  of  infestation. 
Although  all  six  photointerpreters  were  experi- 
enced, none  had  prior  experience  interpreting 
color  infrared  photos  of  lodgepole  pine  killed  by 
the  mountain  pine  beetle,  and  only  three  had 
prior  experience  with  panoramic  photography. 

The  imagery  was  viewed  under  9x  to  12x  mag- 
nification6 with  specially  adapted  microfiche 
viewers  (figure  6).  Quick  estimates  of  the  num- 
ber of  faders  were  made  in  each  cell  having 
lodgepole  pine-type  and  entered  on  a standard 
form  (figure  4).  Cells  estimated  as  having  less 
than  10  faders  (0.06  faders/ acre)  were  recorded 
as  zero.  A tabular  summary  of  the  classification 


phase  follows: 

No.  frames  examined  79 

No.  160-acre  cells/frame  (4  x 38)  152 

No.  acres/frame  (152  x 160)  24,320 

Total  acres  non-type  (3332  cells  x 160)  533,120 

End  lap  (349  cells  x 160)  55,840 

Zero  count  (5073  cells  x 160)  811,680 

>10  faders  (3254  cells  x 160)  520,640 

Total  acres  examined  1,921,280 

(12,008  cells  x 160) 


5.  Scanning  the  OB  imagery  to  eliminate  the  obvious  non- 
countable  cells. 


6 Using  a 9x  magnification,  a single  160-acre  cell 
(at  nadir)  would  fill  the  entire  screen.  With 
increased  magnification,  only  a portion  of  the  cell 
would  be  viewed  at  a time. 
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b-  Form  used  to  record  fader  estimates  during  classifica- 
tion phase. 
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6.  Photointerpretation  team  classifying  the  OB  photo- 
graphy by  1 60-acre  cells. 


Once  this  phase  was  completed,  the  data  was 
arranged  by  computer  program  for  photo  plot 
selection  by  PPS.  In  addition,  individual  cell 
quick  estimates  were  arbitrarily  assigned  dam- 
age intensity  classes  ranging  from  1 to  9 for 
display  in  a map-like  computer  printout 
(figure  7). 

Photo  Interpretation 

The  120,  160-acre  photo  plots,  selected  by 
PPS,  were  examined  with  Bausch  and  Lomb 
240  zoom  stereomicroscopes  (7x  to  lOx)  by 
three  experienced  photo  interpreters  (figure 
8).  Each  photo  plot  was  divided  into  sixteen  10- 
acre  subplots  by  a specially  produced  narrow 
line  (1  mil)  grid  attached  with  transparent 
tape  to  the  back  of  the  frame  (figure  9)7.  Once 
attached,  the  equal-area  grids  on  the  emulsion 
side  were  removed. 

All  faders  within  each  10-acre  plot  were 
counted  and  recorded  (figure  10).  Once  this  stage 
was  complete,  thirty  10-acre  subplots  were  selec- 
ted by  PPS  for  ground  truth.  However,  only  20 
plots  were  actually  visited  and  cruised  because 
of  personnel  and  time  constraints. 

Ground  Truth 

During  past  photo  surveys,  ground  truth  was 
obtained  in  late  summer  or  early  fall  following 
aerial  photography.  Since  film  was  not  pro- 
cessed and  duplicated  until  September  12,  there 
was  insufficient  time  left  to  accomplish  the  pre- 
liminary stages  and  obtain  ground  truth  before 
the  onset  of  winter.  By  spring,  crowns  of  the 
beetle-infested  trees  possibly  could  deteriorate  to 
an  extent  precluding  differentiation  from  the 


8.  Photointerpreter  counting  faders  with  the  stereomicro- 


scope. 

previous  year’s  kill.  This  concern  was  unfounded, 
as  the  faders  retained  their  foliage  through  the 
winter  and  into  early  summer  (figure  11). 

The  10-acre-square  plots  were  located,  deline- 
ated, and  cruised  with  the  aid  of  a specially  con- 
structed stereo  viewer  (figures  12, 13).  Because 
of  personnel  and  logistical  constraints,  these 
steps  were  done  at  two  different  times.  A two-  or 
three-person  crew8  would  first  locate  the  plot, 
and  using  the  stereoviewer  would  locate  and 
mark  the  plot  corners  and  boundaries  with 
colored  plastic  flagging.  A detailed  map  giving 
directions  into  the  plot  was  then  drawn  on  the 
back  of  the  recording  form  (figure  14)  to  aid  the 
follow-up  crew. 

Later  a three-person  follow-up  crew  consist- 
ing of  two  cruisers  and  a recorder  would  string 
the  boundaries  and  run  parallel  string  lines 
within  the  plot  at  varying  intervals  for  orienta- 
tion. All  faders  five  inches  dbh  and  larger  were 
blazed9  and  recorded.  In  some  areas  the  plastic 
flagging  used  to  designate  corners,  boundaries, 
and  access  to  the  plot  was  removed  if  there  was  a 
possibility  that  it  might  conflict  or  be  confused 
with  flagging  established  for  other  purposes  by 
Forest  personnel. 


8 At  times  the  number  of  persons  per  crew  varied, 
but  it  was  felt  that  three-person  crews  were  the 
most  efficient. 


7 Produced  by  the  Forest  Service,  WO-Geometron- 

ics,  Rosslyn,  VA. 


9 The  trees  were  blazed  to  eliminate  double 
counting. 
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7.  Computer  printout  showing  quick  cell  classification  for 
19  alternate  OB  frames.  The  classification  was  arbi- 
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9.  Portion  of  OB  f rame  near  nadir  overlaid  with  1 60-acre 
grid  divided  into  sixteen  10-acre-square  plots  (35mm 
Kodachrome  enlargement). 


10.  Form  used  to  record  fader s by  10-acre  plot. 


11.  More  or  less  typical  lodgepole  fader s on  ground  truth 
plot  No.  109-2-33(B3).  These  trees  were  attacked  and 
killed  by  the  mountain  pine  beetle  in  1977  and  attained 
optimum  fading  in  July  1978.  Photo  (35mm  Koda- 
chrome) taken  May  27. 1979. 


12.  Portable  stereoscope  used  for  viewing  OB  transparen- 
cies on  the  ground  plots.  Most  of  the  detailed  work  was 
done  with  a specially  developed  flat  field,  color-corrected 
7x  stereoscope. 
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Plot  No.  //r-a-f' 

Beaverhead-Gallatin  Optical  Bar  Photo  Interpretation 
Record  for  160-acre  Plot 

Frame  US'  Row  2.  Column  V Interpreter  

Date  Magnification  # 
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13.  Drawing  showing  details  of  portable  stereoscope. 
Inset:  7x  stereoscope  with  detachable  base. 

U.  Ground  truth  tally  form. 
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Estimation  Procedures 


M 

m 
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Stage 
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Variable 

Data  Analysis 

Procedures  for  estimating  the  number  of 
faders  and  total  volume  combine  information 
from  each  step  in  the  design.  These  steps  en- 
tailed (1)  quick  cell  classification  and  sample 
selection  of  cells,  (2)  careful  photo  interpretation 
of  selected  cells  and  sample  selection  of  ground 
plots,  and  (3)  ground  truth  collection.  Sample 
selection  at  each  level  was  based  on  PPS. 


Cell  Classification  and.  Sample  Selection 

Cell  classification,  based  on  the  quick  esti- 
mates, resulted  in  3254  160-acre  cells  contain- 
ing 10  or  more  faders,  or  a total  of  378,339  fa- 
ders. The  first  level  sample  was  selected  by  a 
systematic  PPS  procedure  using  the  following 
parameters: 

Sample  size  = 120 
Total  faders  estimated  = 378,339 
Sampling  interval  = 3153 
Random  start  = 1538 

The  sampling  was  done  with  replacement. 

Carefid  Photo  Interpretation  and  Sample 
Selection 

The  selected  sample  again  was  photointerpre- 
ted.  Data  was  recorded  by  10-acre  subcells.  The 
next  sample  was  determined  based  on  the  ratio 
of  the  careful  count  to  the  quick  count.  A system- 
atic PPS  sample  was  done  based  on  the  following 
parameter: 

Sample  size  = 30 

120 

Total  ratio  = Z Careful/Quick  = 265.61 
i = 1 

Sampling  interval  = 8.85 
Random  start  = 3.55 

The  30  cells  selected  were  done  with  replace- 
ment. One  10-acre  subcell  was  selected  PPS 
based  on  the  number  of  faders  counted  in  each 
subcell. 


Ground  Truth  Collecticm 

The  original  intent  was  to  ground  sample 
thirty  10-acre  subcells.  Personnel  restrictions 
and  time  constraints  limited  the  actual  ground 
samples  to  20.  The  original  sample  was  used  to 
select  20  from  the  30  subcells  selected. 

Data  was  analyzed  on  the  assumption  that 
there  was  a sample  of  20  subcells. 


where:  i = 160-acre  cell 

j = 10-acre  subcell 
k = ground  sample  size 
Qi  = rapid  photo  interpretation  for  cell  i; 
i = 1,  12,008 

M is  the  total  of  all  faders  counted 
Q = Z Qi  from  the  quick  interpretation 
i = 1 

Cjj=  the  careful  interpretation  of  the  jth 
subcell  (10-acre  cell)  from  the  ith 
selected  subcell  (160-acre  cell), 
i = 1,  m;  j = 1, 16 

16  is  the  total  of  all  faders  from 
C;=Z  Cjj  careful  interpretation  from 
i = 1 the  ith  selected  cell. 

Xjj=  ground  measurement  variable 
(trees  or  volume)  from  the  jth 
selected  subcell  and  ith  cell. 

M = total  number  of  cells 
= 12,008 

m = number  of  first  stage  cells  selected 
= 120 

Q =378,339 

n = number  of  ground  subcells  selected 
= 20 


This  formula  shows  the  three  levels  of  sam- 
pling. The  variables  with  single  subscripts  (i),  Qj 
and  Cj  take  on  the  cell  value  for  a given  subcell 
(j).  Similarly,  Xy  is  the  number  of  trees  (or 
volume)  in  the  ground  subcell  (j)  from  selected 
cell  (i). 
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For  computational  ease,  the  formula  can  be 
reduced  to: 

m 

let  Q = £Q;, 

i = 1 


m 

Cj 

R = 

£ 

— then 

i = 1 

Qi 

ii 

[OR] 

f £ ^1/n 

L m J 

L k = 1 Cy  . 

The  corresponding  standard  error  (SE)  is: 


RSE  = — 100,  expressed  as  a percent 

A 

Y 

The  steps  taken  to  derive  individual  expansion 
factors  at  each  stage  in  the  computation  of  plot 
faders  and  volumes  are  shown  in  table  1. 


RESULTS 

Total  lodgepole  pine  mortality  during  1977 
over  that  portion  of  the  outbreak  covered  by  the 
grids  on  the  79  panoramic  frames  (520,640 
acres)  was  estimated  at  1,891,510  trees  and 
27,001,000  cubic  feet  of  volume  (table  2).  Stan- 
dard errors  were  10.3  and  13.6  percent  respec- 
tively. The  multistage  survey  (Bennett  and  Bous- 
field  1980)  covered  270,255  acres  and  reported 
1 ,270,617  trees  and  24,237,335  cubic  feet  of  lodge- 
pole  pine  killed  during  1977  by  the  mountain 
pine  beetle.  These  data,  however,  are  not  compa- 
rable because  of  the  difference  in  area  covered 
by  each  survey.  The  OB  survey  covered  almost 
twice  the  area  (520,640  vs.  270,255  acres)  but  the 
“conventional”  survey  reported  a higher  volume 
per  acre,  indicating  that  the  OB  survey  included 
considerably  more  “light”  infestation.  The  mag- 
nitude of  losses  from  the  OB  survey  would  have 
been  greater  had  not  the  voids  between  succes- 
sive frames  and  sidelap  occurred  (figure  1). 

DISCUSSION 

A regression  analysis,  based  on  the  relation- 
ship between  photo  and  ground  counts  used  in 
previous  double  sampling  and  multistage  sur- 
veys (Wear  et  al.  1966,  Klein  et  al.  1979,  Hos- 
tetler and  Young  1979a),  would  appear  to  be 
inappropriate  for  this  type  of  panoramic  pho- 
tography. The  first  level  of  a regression  analy- 
sis, the  coefficient  of  determination  (R2),  is 
used  as  a measure  of  consistency  rather  than 
accuracy,  and  its  use  is  based  on  the  assump- 
tion that  every  tree  has  an  equal  probablity  of 
being  detected  and  counted.  The  probability  of 
detection  is  unequal,  since  each  panoramic 
frame  has  variable  scale.  For  example,  a tree 
at  nadir  (-scale  1:30,000)  has  a greater  likeli- 
hood of  being  detected  than  a similar  size  tree 
at  38°  (-scale  1:38,000).  Heller  and  Wear 
(1969)  found  that  ground  photo  accuracy  of 
individual  trees  increased  with  scale,  and  that 
a scale  of  1:31,680  was  the  threshold  for  detec- 
tion of  infestations  (groups  of  4-10  trees  per 
group)  smaller  than  50  feet  in  diameter.  This 
probability  of  detection  is  even  further  re- 
duced by  the  obliquity  and  apparent  color 
change  in  the  film  at  increasing  distances  from 
nadir.  These  phenomena  were  first  suspected 


15 


1 (2)  Total  faders  counted  in  Quick  Count  (Q)=  378,339 

(3)  FGrst  stage  sample  size  (m)  = 120 

(4)  Sum  of  ratio  of  careful  counts  to  Quick  count  (R)  = 265.61 
Ground  sample  size  (n)  = 20 
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during  the  photointerpretation  stage  and  later 
empirically  confirmed  during  the  ground 
checks  (all  three  photointerpreters  partici- 
pated in  the  ground  checks). 

An  analysis  was  made  to  determine  whether 
there  was  a change  in  detection  accuracy  with 
distances  from  nadir.  Since  the  survey  was  not 
designed  with  this  objective  in  mind,  the  anal- 
ysis was  more  after  the  fact.  The  twenty  10- 
acre  ground  truth  plots  were  separated  into 
two  10-plot  groups,  the  first  10  falling  between 
nadir  and  ± 19°  and  the  second  10  falling 
between  ± 20  and  40°  (figure  15).  The  ratio  of 
the  ground  count  to  the  photo  count  was  com- 
puted for  each  plot,  arranged  into  its  respec- 
tive group,  summed,  and  a mean  difference 
computed.  The  mean  difference  of  the  ground 
/photo  ratios  were  then  subjected  to  a “T”  test. 
The  difference,  3.67  with  18  d.f.,  was  highly  sig- 
nificant at  the  99  percent  level,  indicating  that 
the  mean  ratios  were  from  two  different  popula- 
tions. Lacking  a more  sensitive  analysis,  it  can 
be  assumed  that  the  photo  counts  were  more 
accurate  within  ± 20°  than  those  beyond  ± 20°. 

Two  important  characteristics  intrinsic  to 
OB  photography,  in  addition  to  variable  scale, 
are  obliquity  and  color. 

Obliquity  may  be  as  important  as  variable 
scale  in  affecting  interpretation  accuracy  and 
the  final  survey  results.  The  angle  that  an  indi- 
vidual vertically  standing  tree  is  imaged  on  the 
film  is  equal  to  the  scan  angle,  i.e.  a tree 
imaged  at  the  30°  point  on  the  film  is  tilted  30° 
from  vertical — away  from  nadir.  In  this  situa- 
tion, faders  in  a closed  or  multistoried  stand 
and  those  on  an  adverse  slope  would  be  either 
totally  or  wholly  masked  by  their  neighbors. 
Small  openings  and  occasional  breaks  in  the 
stand  used  by  ground  crews  for  orientation 
cannot  be  seen.  Codominant,  intermediate,  and 
smaller  trees  become  invisible  in  some  situa- 
tions. 


A significant  color  shift  from  nadir  outward 
occurs.  Faded  trees  at  or  near  nadir  generally 
appear  yellow  to  yellow-orange,  but  become 
less  vivid,  and  in  some  instances  change  to 
tones  of  grey,  with  increasing  angle.  Mountain 
pine  beetle-faded  lodgepole  pine  possibly  could 
have  become  confused  with  Douglas-fir  defol- 
iated by  the  western  spruce  bud wor m (Choris- 
toneura  occidentalis  Freeman)  in  some  areas. 
However,  this  color  attenuation  did  not  occur 
in  a few  locations  where  faded  lodgepoles  were 
open  grown  or  on  the  nadir  side  of  stands  bor- 
dering open  areas.  Within  the  stands,  however, 
faded  trees  appeared  considerably  less  dis- 
tinct. This  phenomenon  is  probably  due  to  the 
relative  position  of  the  tree  in  relation  to  its 
immediate  surroundings,  such  as  a highly 
reflective  meadow  versus  a closed  stand,  and  to 
its  orientation  and  distance  from  the  camera  in 
relationship  to  the  sun. 

Data  gaps  created  within  and  between  flight 
lines  and  topographic  distortion  prevented 
production  of  an  accurate  infestation  map.  The 
flight  was  planned  so  that  successive  frames  on 
parallel  flight  lines  would  overlap  at  the  40° 
point10,  or  approximately  8.5  nm  (9.8  miles) 
from  nadir.  However,  the  aircraft  strayed  off 
course  on  the  center  flight  line  (figure  1 ) result- 
ing in  more  than  sufficient  sidelap  to  the  west, 
but  a significant  holiday  to  the  east.  Unfortu- 
nately, some  of  this  omitted  area  contained 
heavy  beetle  damage.  Holidays,  resulting  from 
insufficient  (1)  overlap  and  (2)  intrack  cover- 
age by  the  equal-area  grids  also  occurred 
between  alternate  frames.  The  overlap  at 
nadir  was  50  percent,  which  also  occurred  in  a 
previous  study  (Klein  et  at.  1979).  The  master 
equal-area  grid  covered  only  two  miles  of  the 
approximate  2.6-mile  intrack  distance  with 

10  This  assumed  mean  datum  at  5000  feet  and  air- 
craft at  65,000  feet. 


Table  2.  Estimates  of  numbers  of  trees  and  volume  of  lodgepole  pine  killed  by  the  mountain 
pine  beetle  in  the  Beaverhead  and  Gallatin  National  Forests,  1977. 


Number  of  Trees 

Percent 

Volume  (M  cu.  ft.) 

Percent 

Acres1 

total  S.E. 

S.E. 

per  tree2  total  S.E. 

S.E. 

520,640 

1,891,510  194,804 

10.3 

0.0143  27,001  3682 

13.6 

1.  Includes  all  160-acre  cells  which  contained  10  or 

2.  Derived  by  dividing  Volume/Trees 

more  faders 

17 


15.  Locations  of  the  twenty  10-acre  ground  plots  and  their 
respective  160-acre  photo  plots  in  reference  to  their 
locations  on  an  optical  bar  frame.  The  20  ground  plots 
fell  on  12  different  frames. 
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insufficient  overlap  preventing  complete  ste- 
reo coverage  of  some  of  the  photo  plots.  The 
grid  coverage  problem  could  have  been  par- 
tially resolved  by  attaching  grids  to  successive 
frames  rather  than  to  alternate  ones,  but  due  to 
the  intrinsic  elevational  distortions,  they  still 
would  not  have  matched.  An  attempt  to  match 
or  separate  individual  160-acre  cells  in  this 
region  would  do  little  to  improve  the  overall 
accuracy.  The  same  problem  would  be  encoun- 
tered in  the  endlap  regions  where  some  cells 
would  be  partially  or  totally  duplicated. 
Errors  resulting  from  aircraft  attitude,  varia- 
tions in  elevation,  and  frame  alignment  com- 
pounded these  alignment  problems.  Some 
frames,  particularly  those  in  the  north  half  of 
the  middle  flight  line,  were  skewed  or  not 
aligned  perpendicular  to  the  line  of  flight.  This 
probably  resulted  from  crosswinds  (crab)  or 
improper  ground  speed11. 

The  equal-area  grid  is  an  extremely  useful 
aid  for  the  transfer  of  information  from  photo 
to  map  and  for  establishing  a systematic  sam- 
pling system.  In  flat  or  even  moderately  rolling 
terrain  some  position  accuracy  can  be  sacri- 
ficed, but  in  mountainous  terrain  where  eleva- 
tions within  a particular  frame  could  vary  as 
much  as  5000  feet  and  where  full  coverage  and 
a relative  accurate  data  base  is  needed,  map- 
registered  grids  should  be  used  when  possible. 
These  grids  compensate  for  aircraft  altitude 
and  displacement  errors  induced  by  elevation. 
They  generally  are  accurate  to  within  ± 200 
feet  but  are  considerably  more  expensive  to 
produce  than  the  equal-area  grids  ($47  vs. 
$1.50).  However,  it  is  believed  this  cost  can  be 
reduced  by  approximately  50  percent  with 
only  a slight  loss  in  accuracy  by  using  photo- 
grammetric  bridging.12. 

The  ground  phase  of  the  survey  went  sur- 
prisingly well.  Most  of  the  difficulty  in  locating 
plot  corners  and  boundaries  were  in  homogen- 
ous one-story  stands  that  occurred  near  the 
outer  edges  of  the  panoramic  frames.  Of  the 
corners  sought,  85  percent  were  found,  8 per- 


11 The  KA80A  OB  camera  is  calibrated  to  operated 
at  constant  speed,  synchronized  with  a ground 
speed  of  400  nm  (460  miles)  per  hour.  Both  sys- 
tems operate  independently.  If  aircraft  speed  is 
less  or  greater  than  400  nm/hour,  skewness  will 
result. 

12  Liston,  op  cit. 


cent  were  judged  as  “possibly”  found  and  7 
percent  were  not  found.  These  corners  were 
located  by  intersection  of  boundary  lines.  The 
stereoviewer  worked  under  all  lighting  condi- 
tions, even  in  the  rain. 

The  cost  of  the  OB  flight  was  $25,000,  but 
this  also  included  panoramic  photography  of 
other  mountain  pine  beetle  infestations  in  por- 
tions of  the  Lewis  and  Clark  and  Flathead 
National  Forests  and  Glacier  National  Park. 
The  cost  of  the  U-2  aircraft  fitted  with  the 
KA80A  camera  was  $3750  per  flight  hour13.  If 
the  Beaverhead-Gallatin  survey  had  been 
flown  as  a single  mission,  the  cost  of  the  photo- 
graphy was  estimated  to  have  been  $16,376. 
This  included  1080  feet  of  film,  one  duplicate, 
processing,  ferry  time  to  and  from  the  project 
site,  and  photography  time  over  the  area.  The 
roundtrip  distance  from  Moffett  Field,  Cali- 
fornia, to  Bozeman,  Montana,  is  approxi- 
mately 1680  miles. 


13  R.  Ekstrand,  Airborne  Missions  and  Applica- 
tions Division,  NASA  Ames  Research  Center, 
Moffett  Field,  CA.  Personal  communication. 
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CONCLUSION  RECOMMENDATIONS 


Number  of  lodgepole  pines  and  volume 
killed  by  the  mountain  pine  beetle  can  be  effec- 
tively measured  by  use  of  OB  photography. 
Past  methods  have  required  two  separate  and 
relatively  time  consuming  steps,  sketchmap- 
ping14 and  aerial  photography.  In  effect,  both 
of  these  information  levels  were  obtained  in 
less  than  an  hour  by  a U-2  aircraft.  The  initial 
classification  phase,  comparable  to  stratifica- 
tion by  sketchmapping,  took  4-6  persons  two 
weeks  to  complete,  while  the  actual  sketchmap 
flight  was  completed  in  less  than  four  days. 
Other  times  are  not  readily  comparable, 
because  the  large-scale  photography  was  de- 
layed until  October  1978  by  contract  and 
weather  problems.  The  OB  survey  would  have 
been  more  efficient  and  resulted  in  a more  re- 
liable map  product  had  there  been  more  accu- 
rate aircraft  navigation,  greater  overlap  be- 
tween flight  lines  and  successive  frames,  faster 
film  processing,  and  use  of  a map-registered 
grid  system. 

A major  concern  of  all  prospective  users  of 
NASA-obtained  aerial  reconnaissance  photo- 
graphy is  its  dependability  and  future  availa- 
bility. If  we  and  other  investigators  continue  to 
explore  its  potential  in  forest  insect  and  disease 
detection  and  develop  it  to  a point  where  it 
could  become  operational,  will  it  be  available 
when  and  where  we  need  it?  Also,  if  the  use  of 
OB  and  other  types  of  high  altitude  reconnais- 
sance photography  by  the  Forest  Service  and 
other  natural  resource  agencies  increases,  will 
NASA  be  able  to  meet  the  demand?  Problems 
relating  to  excessive  turn  around  time  (this 
survey)  and  failure  to  complete  the  photog- 
raphy of  a bark  beetle  outbreak  in  Colorado 
within  the  optimum  time  span  have  already 
occurred  (Dillman  et  al.  1980).  At  this  time, 
private  aircraft  with  high  altitude  capabilities 
are  not  available  and  probably  won’t  be  for  at 
least  another  decade. 


14  Sketchmapping  is  a routine  function  and  still 
would  be  required  for  other  forest  pests  not  vis- 
ible on  the  OB  photography. 


In  the  planning  and  execution  of  future  OB 
surveys  for  measuring  annual  lodgepole  pine 
mortality  caused  by  the  mountain  pine  beetle, 
the  following  additions  or  changes  are  sug- 
gested: 

1.  The  effective  working  angle  of  an  OB  frame 
should  not  exceed  ± 40°  from  nadir.  The 
findings  of  this  survey  indicate  that  tree 
count  accuracy  improves  as  one  approaches 
nadir,  but  the  maximum  angle  at  which 
acceptable  accuracy  occurs  is  unknown. 
Future  surveys  for  lodgepole  pine  should  be 
designed  to  establish  this  threshold. 

2.  Flight  lines  should  be  no  further  than  15 
nm  (17  miles)  apart.  This  will  permit  some 
deviation  from  the  planned  flight  line  and 
still  permit  sufficient  endlap  at  the  ± 40° 
point. 

3.  Photographic  requirements  to  NASA  should 
specify  a minuimum  of  55  percent  forward 
lap  (nadir)  at  mean  datum.  In  this  survey 
portions  of  several  photo  plots,  particularly 
those  that  fell  near  nadir,  were  not  in  stereo. 
Fortunately,  none  of  prospective  ground 
plots  that  fell  in  these  areas  were  selected 
for  ground  truth! 

4.  If  a relatively  accurate  map  product  is 
desired,  photo  to  map  transfers  should  be 
made  with  map-registered  grids.  A photog- 
graphic  bridging  technique  should  be  at- 
tempted, although  some  accuracy  may  be 
sacrified. 

5.  Inexpensive  monoscopic  projection  equip- 
ment for  simultaneous  viewing  of  larger 
areas  (at  least  640  acres)  with  at  least  8x 
magnification  should  be  developed. 

6.  The  ideal  crew  size  for  ground  truth  data 
collection  is  three,  with  one  person  record- 
ing and  the  other  two  detecting  and  measur- 
ing the  affected  trees.  The  use  of  visible 
string  lines  is  mandatory  for  definitive 
boundary  separation  and  crew  orientation. 

7.  The  turnaround  time  between  photogra- 
phy, processing,  and  delivery  to  the  user 
should  not  exceed  10  days. 

8.  The  10-acre-square  ground  plot  is  satisfac- 
tory. The  size  and  the  number  of  trees  found 
on  a plot  does  not  seem  to  be  as  time- 
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consuming  as  originally  thought.  More  time 
is  spent  traveling  to  the  plot,  locating  it,  and 
establishing  its  boundaries  than  cruising. 
Of  the  total  plot  time,  52  percent  was  spent 
on  corner  and  boundary  location  while  48 
percent  was  spent  cruising. 
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